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Abstract. - A Landau-Zener multi-crossing method has been used to investigate the tunnel 
splittings in high quality Mni2-acetate single crystals in the pure quantum relaxation regime 
and for fields applied parallel to the magnetic easy axis. With this method several individual 
tunneling resonances have been studied over a broad range of time scales. The relaxation is 
found to be non-exponential and a distribution of tunnel splittings is inferred from the data. 
The distributions suggest that the inhomogeneity in the tunneling rates is due to disorder that 
produces a non-zero mean value of the average transverse anisotropy, such as in a solvent disor- 
der model. Further, the effect of intermolecular dipolar interaction on the magnetic relaxation 
has been studied. 



Single molecule magnets (SMMs) such as Mni2, Fes and Mn4 have enabled new studies 
of the physics of magnetic nanostructures. These materials are ordered crystals of nomi- 
nally identical high spin (S = 10 and smaller) molecules with predominately uniaxial mag- 
netic anisotropy. A particular focus has been on quantum tunneling (QT) of the magnetiza- 
tion |,||,@, |, with efforts aimed at understanding the origin of the tunneling 
phenomena in SMMs, the effect of intermolecular interactions in SMM crystals and the role 
of the environmental degrees of freedom in modulating QT (including phonons and nuclear 
spins) p| . Even in the original, high symmetry and high magnetic anisotropy barrier SMM, 
Mni2-acetate, the QT phenomena remains poorly understood. Primarily this is due to the 
fact that the magnetic properties do not reflect the tetragonal symmetry (S4 site symmetry) 
of the molecule. For example, QT transitions do not obey the selection rules imposed by this 
symmetry and the magnetic relaxation is clearly non-exponential, indicating inhomogeneity in 
the tunneling rates. It has become evident in recent research that disorder that lowers 
the molecule symmetry plays an essential role in the physics of QT in Mni2-acetate single 
crystals, but the precise nature of the disorder remains unclear. Thus far two alternative 
models of disorder have been proposed. The first is line-dislocations in crystals that produce 
long-range deformations and through magnetoelastic interactions lead to a distribution of the 
transverse magnetic anisotropy. This leads to a very broad distribution in the tunnel splittings 
on a logarithmic scale . The second is disorder in the solvent molecules that surround each 
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Mni2 core that lead to a discrete set of transverse anisotropics |15 . Each model has predic- 
tions that can be compared to experiment. In addition, intermolecular dipolar interactions 
are another important source of magnetic inhomogeneity that can be studied in experiment. 

In this letter we present experimental results on the tunnel splittings in deuterated Mni2- 
acetate single crystals in the pure quantum regime, in which relaxation is by QT without 
thermal activation jl^. Deuterated crystals have been studied because the purity of the 
chemicals used in the synthesis leads to very high quality crystals We have used a 

modification of the Landau-Zener (L-Z) method and a high sensitivity micro-Hall effect mag- 
netometer to study the magnetic relaxation in Mni2. With this method we have determined 
the distribution of tunnel splittings for a single QT resonance in a manner that enables a 
critical comparison of the experimental data with models of disorder for Mni2. We also show 
the effect of intermolecular dipolar interactions on the magnetic relaxation in Mni2. 

The effective spin Hamiltonian for Mni2 is 

n = -DSl-BSt-gfiB'ii-S + H' , (1) 

where the parameters D = 0.548(3)K, and B = 1.17(2)x 10^"^K, have been determined by 
EPR and inelastic neutron spectroscopy experiments j^. The first two terms represent 
the uniaxial anisotropy of the molecules, the third term is the Zeeman interaction of the spin 
with the magnetic field, H, and the fourth term includes dipolar and hyperfine interactions 
and higher order transverse anisotropy terms that do not commute with S^. The observed 
steps in the magnetic hysteresis correspond to an enhancement of the relaxation towards the 
equilibrium magnetization when the eigenstates of 5*2, m and m', on opposite sides of the 
barrier arc nearly degenerate (see, for example. Fig. 1 of Ref. These resonances are 

labeled by the integer k = m ~ m' and occur at longitudinal magnetic fields of approximately 
Hk = kD/{giiB) - k 0.42T. 

It has been shown that the L-Z method is a powerful way to study tunnel splittings in SMM 
crystals The method consists of sweeping the longitudinal magnetic field at a constant 
rate, a = dH/dt, across a QT resonance and measuring the fractional change in magnetization. 
For an ensemble of non-interacting identical SMMs, this fractional change in magnetization 
is related to the probability that an individual SMM has remained in the original metastable 
state after an avoided level crossing. The probability is R = {M — Meq) / {Mini — Meg), where 
Mini is the magnetization before the crossing, Al^q is the equilibrium magnetization, and M 
is the magnetization measured after crossing the resonance. This probability is related to 
the quantum splitting A, through the L-Z formula Riz = exp(— 7rA^/2i/oa) , where i^o — 
gfisi^S — k) and I'oa is the energy sweep rate. Clearly, for this formula to be valid the energy 
sweep rate that an individual SMM in a crystal experiences must be proportional to the sweep 
rate of the applied external field. This will only be true if changes in the internal magnetic field 
on each crossing are small {R w 1). For this reason, in our experiments, we cross a resonance 
(say, the fc-th resonance) n-times at a constant and sufficiently fast sweep rate, so that on 
each crossing R k \. Then we repeat this process for different sweep rates. In this case, the 
probability to remain in the metastable state after the n-th crossing is Rk = , where 

Mn is the magnetization after the n-th crossing of the resonance. The L-Z probability after 
n-crossings is: 

Rizn = exp ~- , (2) 

V 2(^0 aj 

Based on this formula, relaxation curves recorded at different sweeping rates should scale when 
plotted as a function of a^ff — a/n. In figure 1 we show the data obtained for resonance 
A: = 7 at different sweeping rates, a = 3.33x10"^, 6.66x10"^ and 13.3xlO~^T/s in the pure 
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Fig. 1 - The tunneling probability for deuterated Mnl2-acetate as a function of the effective sweeping 
rate, a^ff = a/n of the applied magnetic field in a L-Z multi-crossing experiment recorded at different 
sweeping rates (small closed symbols) for resonance k =7 and at a temperature of 0.65 K. The open big 
squares correspond to the tunnel probability in a L-Z single-crossing experiment versus the sweeping 
rate (n =1). The difference between the two curves we associate with dipolar interaction between 
molecules llStl. 



quantum tunneling regime for this resonance, T = 0.65K (to = 10, to' — -3). To do the 
measurements we saturated the system with a negative high magnetic field, H = -5T, and 
then we swept the field to positive values until we arrived to the vicinity of this resonance 
{Hr = 3.55T), with the applied field aligned with the z-axis of the crystal . Then we crossed 
the resonance n times ramping the field up and down within ±0.2 Tesla of the resonance at 
a given rate. This procedure was then repeated for different sweep rates. The scaling of the 
curves is clear, indicating that ae// is the relevant parameter to characterize the data for this 
range of sweep rates. It is important to point out that the relaxation curve does not follow 
eq. (^, i.e., it is not an exponential function of ae//, which we discuss in detail below. 

For smaller sweep rates, changing internal magnetic fields should affect the magnetic relax- 
ation. This was studied experimentally in Fes [@ , where the relaxation departs from the L-Z 
behavior for sweep rates below ~10~^T/s, and was theoretically explained in terms of dipolar 
interactions by Liu et al. To study this in Mni2, we have measured a single-crossing L-Z 
relaxation at resonance fc = 7 for different sweeping rates ranging from ~10~^ to ~10~^T/s. 
The results are compared to the multi-crossing procedure in fig. 1 (open squares). The differ- 
ence between these procedures is clear. The break point from the L-Z single-crossing behavior 
is ckc ^ 10~^T/s and is similar to that observed in Feg, while the deviations from the scaled 
data are smaller. This is due to the fact that dipolar interactions in Mni2 are smaller than 
in Fe8 because of differences in the intermolecular distances and crystal structure We 
emphasize that the L-Z multi-crossing procedure permits studies at very small effective sweep 
rates in a regime in which a single crossing experiment would be complicated by intermolecular 
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Fig. 2 - (a) Probability to remain in the metastable well in a L-Z multi-crossing experiment using 
different sweep rates (3.33xlO"^T/s to 1.33x 10"^T/s) for different resonances, k. Lines are fits using 
a log-normal distribution of splittings, (b) Corresponding splitting distribution functions extracted 
using the fitting procedure described in the text. The centers, Xc, are -7.34, -6.54, and -6.07 and the 
widths, W , are 0.43, 0.36, and 0.24, for resonances, k =6, 7, and 8, respectively. 

dipolar interactions. 

We have carried out multi-crossings measurements of resonances fc = 6, 7, and 8 at different 
sweep rates, a > ac- The results are presented in fig. 2a starting from a saturated sample. 
Small differences were observed in the same procedure carried out on three different single 
crystals synthesized in the same way. The main results are the following: (a) The behavior 
of the relaxation is clearly non-exponential. The relaxation is broad on a logarithm scale, 
indicating the presence of a distribution of quantum splittings in the crystal, and, (b) the 
data give direct information on the distribution width due to the fact that multi-crossings 
procedure allows us to measure a large fraction of the relaxation curve for each resonance. 
We have assumed a log-normal function to extract the distribution of the quantum splittings 
in Mni2. We take the form, f{x) — Aexp{—{x — Xc)/2W'^), with x = log A, and we fit the 
relaxation curves with the following expression, 

Ri,{a,A,n)f{logA)dlogA , (3) 

-oo 

The fitting parameters are the center of the distribution Xc and the width of the distribution 
on a log scale, W. The fits with eq. are represented by lines in fig. 1 and fig. 2a and arc in 
excellent accord with the experimental data. The fitting parameters can be determined within 
an accuracy of 5%. The corresponding splitting distribution functions for each resonance are 
shown in fig. 2b. The mean value of the distribution increases with k while the width remains 
more or less constant, being somewhat narrower for k = 8. This width is a factor of 2 (in log A) 
smaller than the width extracted in ref. |l^, because of the higher quality of the deuterated 
crystals. 
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Fig. 3 - Splitting distributions function extracted using (A) a log-normal distribution, (B) the disloca- 
tion model, and (C) the solvent disorder model. Distribution (D) was measured in a non-Deuterated 
Mni2-acetate single crystal (from Ref. The inset shows the distribution of the second order 

anisotropy parameter, E, extracted from the log-normal distribution of tunneling splittings using 
eq. for k — 6 and k — 8. 



Recently it has been proposed that line dislocations or solvent disorder [T^] should 
introduce a second order anisotropy term, E{Sl — Sy) in the Hamiltonian of Mni2, that 
increases the quantum tunneling probability. Line dislocations generate a broad distribution 
of the parameter E with mean value < E > = 0, while disorder of the solvent molecules creates 
a discrete number of E values due to different Mni2 isomers existing in the crystal. Recent 
relaxation experiments have confirmed the existence of a broad splitting distribution [p2|,p^ 
as well as broad absorption lines ||2l| and discrete double absorption peaks jlj, ^have been 
observed by EPR experiments. In order to compare our experimental data to these models 
we have tried to fit our relaxation curves to the distributions predicted by these models. In 
the distribution function predicted by the line dislocations model the mean value and width 
of the distribution are not independent variables, 

where x = hiE with E = E/2D. Eg is the width of the distribution of the anisotropy param- 
eter E. Ec depends on the geometry of the crystal and on the concentration of dislocations 
per unit cell, c. To fit to this function we use the following relation between A with E which 
follows from perturbation theory, 

ln(Afe/.g,)/6=ln(^^) , (5) 

where gk and £,k are determined by k, S, and D mm. The only fitting parameter is Ec- The 
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result of this fit for resonance A: = 6 is shown in fig. 3 (dashed-dot hne). It corresponds to a 
concentration of dislocations of around c 10^*. The mean value is chosen to coincide with 
the center of the log- normal distribution function (continuous thick line) . Clearly, the width is 
many orders of magnitude too large to fit the experimental data. Physically this is due to the 
fact that the transverse anisotropy distribution expected from the dislocations model alone is 
centered at < > = 0, which produces a large tail to small values of the tunnel splitting, 
and would produce a much broader relaxation curve than observed in experiment. As a 
consequence it is clear that a model that describes our data must have non-zero mean value 
of the transverse anisotropy distribution. For example, taking our log-normal distribution of 
A, the distribution of the transverse anisotropy parameter can be inferred. The results for k 
= 6 and k = 8 are shown in the inset of figure 3. As these curves do not overlap (scale) it is 
clear that the origin of the tunnel splittings cannot only be due to a second-order transverse 
anisotropy. However, qualitatively, this would lead to a distribution of A that narrows for 
increasing k (which can be seen from perturbation theory, eq. (||)), as observed in fig. 2. 

The solvent disorder model suggests a multi-peak discrete splitting distribution due to 
the presence of six different isomers of Mni2 with different E parameters. We have found 
that it is not possible to fit the experimental data with a sequence of delta-functions for the 
suggested values of E for each isomer. However, we have obtained a good fit to the data 
assuming a log-normal distribution with a very narrow width for each delta-function proposed 
in ref . . The resulting distribution function has three well defined peaks as can be observed 
in fig. 3 (continuous thin line). The values of the peak centers are x^i =-7.19(-7.0525), 
Xc,2 =-8.55(-8.1749), and Xc.s =-6.60(-6.7995), in accord with the reported values between 
parenthesis The width of each peak is Wi = Xd/^O- The height of each peak is taken to 
be proportional to the measured population of its corresponding isomer. The data suggest that 
the solvent disorder model, which gives a non zero value of the transverse anisotropy, together 
with some other type of disorder (like line dislocations or point defects) that introduces a 
small broadening of the anisotropy E parameter around these values are consistent with our 
relaxation data. We note other models that include higher order transverse anisotropics 
(such as C(S'^ — Sy)"^) and disorder of some type will likely also have features necessary to 
understanding the relaxation in Mni2. Of course, in order to explain the relaxation for odd 
resonances there must be transverse fields due to small tilts of the anisotropy axis of the 
molecules or slight misalignments of the applied field. 

To conclude, we have used a powerful method to study quantum splittings of SMM single 
crystals. We have reported for the first time the break-down of L-Z scaling due to dipolar 
interactions in Mni2-acetate and we have estimated the order of magnitude of the interaction 
energy between molecules. Our results suggest that the quantum splitting distribution in 
Mni2-acetate should be due to disorder which generates a transverse anisotropy distribution 
with a non-zero mean value, such as solvent disorder. 
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